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Abstract 

We have developed a model for the NN —> NNinr reaction and evaluated 
cross sections for the different charged channels. The low energy part of those 
channels where the pions can be in an isospin zero state is dominated by N* 
excitation, driven by an isoscalar source recently found experimentally, followed 
by the decay N* —> N(irir)'g= ( ) )ave . At higher energies, and in channels where 
the pions are not in T=0, A excitation mechanisms become relevant. A rough 
agreement with the experimental data is obtained in most channels. Repercus¬ 
sions of the present findings for the ABC effect and the pp —* ppn 0 reaction close 
to threshold are also suggested. 
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1 Introduction 


Pion production in NN collisions is one of the sources of information on the NN inter¬ 
action and about nucleon resonance properties. Particularly the two pion production 
channel might be especially enlightening in view of the interesting information obtained 
from the study of analogous reactions with two pions in the final state, the tiN —> 7T7 t N 
and yiV —*■ nnN reactions. 

The nN —> nnN reaction close to threshold has been a testing ground for chi¬ 
ral symmetry [jlj in the tiN sector, although as one goes away from threshold the 
contribution of resonances becomes important 0. Even at threshold there is an im¬ 
portant non-vanishing contribution from N* (1440) excitation followed by the decay 
N* —>■ N{'n'n)s- i l a v e [0- With the formalisation of chiral perturbation theory one can 
further exploit the ideas of chiral symmetry and perform calculations including loops 
with effective Lagrangians Q]. The Roper resonance still plays and important role but 
one can use effective Lagrangians involving resonances up to order 0(p 2 ), which gener¬ 
alise the coupling of the N* —> N (nn)'g=^ ave transition used in Ref. |2| and introduce 
one new degree of freedom into the scheme 0 . 

In the ylV —> nnN reaction, which has also been the subject of study at threshold 
as a further test of chiral symmetry [f|, [j|, one of the most interesting results is the 
role of the 1V*(1520, J n = 3/2“) resonance |J. The mechanism of N* (1520) photo¬ 
excitation followed by the decay into Air interferes with the dominant Kroll Ruderman 
A Nn'y term and allows one to obtain information about the N* (1520) —> An decay 
amplitudes, additional to the one obtained from the analysis of the nN —» nnN reaction 
0, and which poses new challenges to quark models of baryons ||. Roper excitation 
also plays an important role at threshold in the 7 IV —> nnN reaction but its role is 
shadowed by the contribution of other resonances at higher energy |J . 

A common feature of the nN nnN and 7 N —> nnN reactions at threshold is 
the role of the N* (1440) followed by its decay into N {nn^z^ave ■ This gives us a 
hint that Roper excitation might also play an important role in the NN —> NNnn 
reaction close to threshold. On the other hand, contrary to the case of N* excitation 
from 7 tN and 7 N , which is well known, now we have to deal with the NN —> NN* 
transition which is not so well explored as the one baryon vertices. An important step 
in understanding the NN —> NN* reaction was given with the experiment of Ref. 
m which showed a relatively large strength for N* excitation due to the exchange of 
an isoscalar source, the only one allowed in the (a, a') reaction on the proton when 
this proton is excited to the Roper. The experiment has been analysed in Ref. Jll 


where the large background from delta excitation in the projectile and an important 
interference term between the latter mechanism and the one from Roper excitation 
in the target are properly accounted for. The analysis provides the strength for the 
isoscalar NN —> NN* transition which turns out to be large compared to ordinary one 
pion exchange. 

The combination of this N* isoscalar excitation with the N* —>■ N (nn^z^ave decay 
seems then called to play an important role in the NN —>■ NN nn reaction and this 
will be the case, as we shall see. On the other hand in the NN —> NN nn reaction 
there are many terms which contribute, starting from the set of terms provided by 
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chiral Lagrangians, plus terms mediated by A resonance excitation. We study all these 
terms asserting their relevance as a function of energy. Furthermore we observe that 
the weight of the different terms varies appreciably from channel to channel and hence 
the combined information from these channels puts strong tests to the model. 

The comparison with the limited experimental information available on this reaction 
0300 shows that a rough agreement in most channels can be obtained. Remaining 
discrepancies in some channels are discussed. 


2 Contribution from non-resonant terms 


The model that we use is depicted diagrammatically in Fig. 1 and contains terms 
proceeding through baryon resonance excitation and other non-resonant terms. The 
latter ones are shown in diagrams (1), (2), (3) and require the knowledge of the ttN —> 
7 tN amplitude, the hNN vertex, the nn —> tttt amplitude and the NNiririr vertex. 
In the spirit of the phenomenological approach followed here, we take the empirical 
7 tN —> 7iN amplitude. The other ingredients are derived from the lowest order SU(2) 
chiral Lagrangians at tree level [|1^, 16[ . 

The lowest chiral Lagrangians which we need are given by 


C = C 2 + C[ B) ( 1 ) 

where £ 2 is the Lagrangian in the meson sector and contains the baryon sector 
with the meson baryon interactions. We have 

f' 2 

Co = < d^d^U + M{U + W) > (2) 

C\[ B) = - M b + |V75^ , (3) 

where the symbol <> denotes the flavour trace of the SU(2 ) matrices and U, u are 
2 x 2 matrices 


= w(0) 2 = exp (iv / 2 < L//) 


where 


(1> ( X ) 


= ( 72 71-0 7r+ 

1 71 ~ 72 n ° 


and 


T = 


P 

n 


The mass matrix M appearing in Eq. (2), in the limit of m u = rrid, is 


(4) 

(5) 

( 6 ) 
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(7) 


M 


ml 0 \ 

0 ml J ' 


Furthermore, V M is a covariant derivative defined as 


= d^i + , 

1 % = \{u^d^u + ud^v)) , 


( 8 ) 


= iv)dJJv) . 

By expanding the Lagrangian £ 2 up to order 0(</> 4 ) and keeping the interaction terms 
in C[ B ^ up to order O(0 3 ) as done in Ref. [T?|, we obtain the diagrams (1), (2) and (3) 
of Fig. 1 which contribute to the NN —> NNnn reaction. 

The term with in C\ ; contributes to the isovector part of the s-wave ttN ampli¬ 
tude. The isoscalar part of this amplitude is zero at this order in the chiral expansion 
and empirically is much smaller than the isovector one. The isoscalar amplitude is 
generated when terms of order 0(p 2 ) in the baryon chiral Lagrangian are introduced 
0. In this latter reference the couplings of the new Lagrangians are deduced from 
empirical information of the nN amplitude in Ref. [|3|. With the same philosophy we 
take the isoscalar nN amplitude from experiment and add it to the vector part. It is 
customary to write the isoscalar nN —> nN amplitude in terms of Ai. Using Mandl 
and Shaw normalisation (R| ( —it = Ad of Mandl and Shaw ) one has 


2A 

— it^ = —iAn —- ; Ai = 0.0075 ; /i = m n (9) 

h 

which is diagonal in spin and isospin. The isovector part for pn + —> pn + amplitude is 
given by 

- it^ = -i-^—pYp^p^ + j/J (10) 

where p,p' are the initial, final pion four-momenta and p (p) besides the Y matrix 
stand for the spinor of the final (initial) protons. This would compare with the usual 
empirical form of the low energy amplitude 


— it^ = —iAn—^ip + p) 0 ( 11 ) 

b 

which comes from Eq. (10) by retaining the dominant zero component, with the 
equivalence A 2 = /i 2 /167r/ 2 , which holds up to 15%. We divert from the standard 
chiral expansion and use the empirical nN amplitudes in our evaluations. 

Eq. (10) is trivially modified in the other isospin channels: same factor in nn~ —> 
nn~ , relative minus sign in 7r + n —> n + n and n~p —> n~p, zero contribution for 
n°p(n) —> n°p(n), extra factor y/2 for n ~p —*• n°n and n°n —> n~p and extra fac¬ 
tor — y/2 in n°p —> n + n and 7 r + n —> n°p. 

The amplitudes nn —> nn come from the expansion of £ 2 up to order O(0 4 ). In 
the diagram (1) of Fig. 1 they appear with two pions off shell. We need four charged 
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amplitudes from where we deduce the rest by using crossing symmetry. We shall give 
the explicit expressions only in the 7r°7r° —> 7r + 7r - reaction. The expressions for the 
other channels can be obtained in a similar way. We get 


— it = i 


1 

W 2 


(3 ml + 4 p + ■ p +2 q ■ q) 


( 12 ) 


where q, q' are the 7r° momenta and p + ,p~ the momenta of i r + and tt~ respectively. 

The other vertex needed is the irNN which comes from the ()a terms of Eq. (3) 
expanding up to 0(0). For ir°pp and an outgoing ir° with momentum q it is given by 


- it = Yj irflJ ^ 5 pqi‘ ( 13 ) 

which in the non-relativistic description is — f nNN ffq, the standard Yukawa coupling, 
with the equivalence < 7 ^/ 2 / = /^nn/p- 

Finally the other ingredient needed is the three pion vertex which comes from the 
expansion of the ()a term in Eq. (3) up to order O(0 3 ). For 7T°p —* pn + n~, with q' the 
7T° momentum we obtain 


~ it = ^p 7 M 7s:p _J_(2p+ + 2 Pfl + AcQ (14) 

and similar expressions in other isospin channels. 

With these ingredients one constructs immediately the contribution from diagrams 
(l)-(3) of Fig. 1 for the pp —> ppir + 7r~ case by following the standard Feynman rules. 
As an example, the term corresponding to the diagram (1) of Fig. 1 for pp —> ppir + ir¬ 
is given by 


—it = + Ap + ■ p +2 q ■ q') 


5 ?- T 

q z —m%. q ,z —m £ 


'9Al\2 
, 2 fJ 


{pi^lbPqn) i {pi v lbPq' v 


(15) 


v) 2 


where the indices 1 and 2 refer to the first and second nucleons respectively. In the Ap¬ 
pendix we further reduce this expression to a Pauli spinor notation and care explicitly 
about baryon antisymmetry, the symmetry of the diagrams and off shell form factors. 


3 A resonance terms contribution 

Next we introduce terms with A excitation. We need the irNA coupling, whose vertex 
is given in a Pauli spinor notation by 


— ibH-nNA = —S^qT^ x + h.c. (16) 

p 

corresponding to the vertex ir x + N —> A, with A the isospin index of the pion, q its 
momentum and S' ( T' ) the spin ( isospin ) transition operators from 1/2 to 3/2 
with the normalisation 
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(17) 


2’ 


M 


Si 


, m ) = C 


1 3 

2> g ^ M 


where z/ is the spherical component of and same normalisation for T 7 
We shall also need the 7 tAA vertex for nA —> A given by 


-iSH„AA = —S A qT* (18) 

h 

where now S'a ( 7a ) are the ordinary spin ( isospin ) matrices for spin ( isospin ) 3/2. 
Two properties involving sums over intermediate A spins are needed in the evaluation 
of the diagrams and they are 


Si | M s ) (M\ Sj = \Sij - ^e ijk a k 
m 3 3 

Si I M') (M'\ S A ,j | M) (M I sl = ^ie ijk - ^S tJ a k + - \bjkVi 

m , m ' b b d b 


(19) 

( 20 ) 


For the coupling constants we take, /* 2 /47 t = 0.36, f A /fnNN = f, where the first 
one is empirical and the second one comes from the quark model. The couplings in 
Eq. (16) (18) have to be evaluated in the rest frame of the created A. 

With these new ingredients and those of the former section one can evaluate the 
diagrams (9)-(15) which appear in Fig. 1. The expressions for the amplitudes of the 
relevant terms are written in the Appendix. As we shall see, all these terms are small 
except diagram (12), the one involving the excitation of a A in each nucleon. This 
diagram is only relevant at energies T p > 1 GeV in the laboratory frame, where the 
terms involving crossed A are small ( diagrams (13)-(15) ). At lower energies all of 
them become negligible since they involve two p-wave couplings of the pions ( diagrams 
(10) and (11) involve only one p-wave coupling but they are small anyway since the 
s-wave coupling is comparatively small ). 

One warning is however in order with respect to diagrams (9), (12)-(15) in Fig. 1 
which involve two p-wave couplings. When this is the case one must take into account 
the indirect effect of the rest of the NN interaction, particularly the repulsive force 
at short distances. This is so because the combination of the pion propagator and the 
p-wave couplings contains a 8(f) force which becomes inoperative in the presence of 
short range correlations. We follow a standard procedure to account for this which is 
to substitute 


mjD(q) -► Vi(q) q4j + VJ.(?)(iy - Mi] 


( 21 ) 


which substitutes the pure spin longitudinal pion exchange by a mixture of spin longi¬ 
tudinal and transverse contributions. Following Ref. |L9| we take 
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Vl(q) = V;(q) + g[(q) , 

Vt(q) = v p(<i) + sfM) > 

VM = F 2 (q 2 )q 2 D 7T (q) , 
V'Aq) = F 2 (q 2 )q 2 D p (q)C p , 


( 22 ) 


sl(?) = k* 2 + k 2 )^A, 


2 2 


3 g 2 F p 2 L) p C p , 


q'M) = -\£h 2 D« ~ (q 2 + lq 2 c )F p 2 D p C p 

where D^^q) are the meson propagators and F ntP (q) are meson form factors of the 
monopole type with = 1.3 GeV, A p = 1.4 GeV and C p = 3.94, consistent with the 
Bonn model |2(J . The corresponding functions with a tilde are obtained by substituting, 
in the argument q of the function, q 2 by q 2 + q 2 , where q c = 780 MeV , the inverse of 
the short range correlations scale. 


4 iV*(1440) terms contribution 

Next we introduce the N* (1440) excitation. The N* couples to Nti , Att , N 
as important decay channels. The corresponding vertices are given by 


/ ^ 

~i5H ntNN * = —a qr x + h.c. 
A 

■ c Tj Q'kN *A n t 

- %oH nN * A = - 5 'q T 

h 


(23) 

(24) 


both evaluated in the N* rest frame. The coupling constants are / = 0.477 and 
9ttN*a = 2.07 0. 

For the N* —> N (7T7r);gC° a? , e decay channel, which accounts for 5 — 10% of the total 
width according to the particle data table |2l|, we take the 0(p 2 ) Lagrangian as given 
in Ref. || using an SU(2) formalism 


F n *Nittt = c i'0JV*X+V’iv - 77 '^{d p d u 'ilj N *)u tl u v ^ N + h.c. 


M 


(25) 


where 


X+ = h 2 ( 2 - + •••); u » = -jrdni + ... (26) 

and (j) stands for the pion held. Hence, expanding up to second order in the pion fields 

2 

and keeping the largest /x, v — 0, 0 components, which neglects terms of order 
we find 


- (j) 2 - 1 

F N *Nwn = + C* 2 'ljj N * — (fd 0 (l))(Td 0 (p)'llj N + h.C. 


f 2 


f 2 


(27) 
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This Lagrangian leads to an effective vertex 


- i5H N . Nmr = -i2^(c{ + cp^\ (28) 

for 7t + ti~ and 7 r° 7 r° creation and zero otherwise. Here are the energies of the 

pions. 

Following the steps of Ref. H we fit the experimental width from 7 r° 7 r° plus 7 r + 7 r - 
final states 


^ — ®(Ci) 2 + P(c^) 2 + Z c l c 2 


(29) 


and we take T jv * Nn7! assuming T* N * tot = 350 MeV |2l| and branching ratio of 7.5% for 
the N(7nr)s-° ave decay channel. We obtain a = 0.49710” 3 GeV 3 , [3 = 3.66 10 ” 3 GeV 3 
and 7 = 2.69 10 “ 3 GeV 3 . These values are somewhat different from those of Ref. |3j 
where a smaller width for the Roper was taken. With the c* parameter alone one gets 
a solution for that parameter, but with the two parameters one gets an ellipse in c*, 
c 2 . This ellipse spans over a large range of values of c*, c 2 . It is clear that one needs 
further constrains in order to determine c ^ and c 2 . We use the tt~ p — > 7r + 7r _ n reaction 
in order to impose such constrains. For this purpose we use the model of Ref. 0. 
The cross section shows a strong sensitivity to the values of c 2 and spans about 
two orders of magnitude when the values of these parameters are varied along the 
ellipse of Eq. (29) ( details can be seen in Ref. m )• The best agreement with the 
experiment is obtained with c* 2 = 0, which corresponds to cl = —7.27 GeV^ 1 ( set I ). 
However, the experimental errors are compatible with the use of c\ = —12.7GeR _1 , 
C 2 = 1.98 GeV^ 1 ( set II ) and intermediate values in the ellipse. I 11 order to illustrate 
the uncertainties in c*, c 2 we plot the cross sections using the extreme values of these 
parameters, compatible with the tt~ p —> TT + TT~n cross sections, quoted above. A best 
fit to the different nN —> ttttN channels is done in |23| using a Lagrangian for the 
N* —> N{TnT^Zyja^, process, equivalent to the one of Eq. (27) at low energies. Although 
the data would roughly be compatible with the use of C 2 = 0 , they obtain a best fit 
with other set of parameters. However, this work contains some small differences with 
respect to the model of Ref. [P|, like the relativistic treatment, the use of an N* width 
assuming phase space for nN decay only and the neglect of the isoscalar nN amplitude 
and intermediate A A states ( which are included in |2| ). These are small differences, 
but they have some influence in the and parameters. For the reason of consistency 
with our input, we take the values of c* and c* 2 from our own analysis. 

The diagrams which we obtain now with these new ingredients are (4)-(8) of Fig. 
1. Diagrams (4) and (5) are relevant but ( 6 ) and (7) give a much smaller contribution. 
Diagram ( 8 ) is also important because both N* and A can be placed simultaneously 
on shell. Note that we disregard terms like ( 8 ) or (9) with crossed N* or A poles. Also 
in the case of the N* we disregard the term with two N* excitations, both in the same 
nucleon or in different ones. Because the nNN coupling is about one fourth of the 
nNA one, and the N* propagator would be placed on shell at higher energies, these 
terms are small compared to the corresponding ones with two A or N* A excitation in 
the energy range where we are. 
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In Fig 1. we have explicitly shown what kind of particle exchange we consider. For 
instance, in diagrams (6), (7) we exchange one of the pions coining from the N* decay. 
In diagrams (4), (5) and (8) we consider the effective isospin T=1 interaction of Eq. 
(21), but in addition we must consider an exchange in the T=0 channel, to which we 
come below. 

The analysis of the (ct, a') reaction on a proton target 


TDI carried out in Ref. 


0 interpreted the results in base of two mechanisms: A excitation in the projectile 
depicted in Fig. 2 (a) and N* (1440) excitation in the proton target, Fig 2 (b). 
However, an important interference was found between the mechanism of A excitation 
in the projectile and N* excitation in the target followed by the decay of the N* 
into NTi. Since the 4 He beam has T = 0 the N* excitation in the (a, a') reaction 
requires the exchange of an isoscalar object. In meson exchange pictures it could be 
cr or to exchange. However, the experiment has not enough information to provide the 
separate strength of both ingredients and only the strength of the combined exchange 
can be extracted. This transition amplitude was parametrised in Ref. [|H[] in terms of 
an effective V” which couples to NN as the a exchange of the Bonn model |20| and 
couples to NN* with an unknown strength which is determined by a best fit to the 
data. Hence, we have 


- i5H aNN = F a (q)g a NN ; gl NN /4n = 5.69 
—idHaNN* = F C7 (q)g a NN * 


(30) 


where F a {q) is a form factor of the monopole type with A a = 1700 MeV, assumed 
equal in both vertices. The fit to the data provides a value g% NN * / 47T = 1.33, with 
9<tNN and g a NN * of the same sign. 

In our selection of the relevant diagrams we have excluded those with a nucleon 
pole. These terms are smaller than those with A pole because the nNN coupling is 
about a factor two smaller than the nN A one, and also because there are usually large 
cancellations between direct and crossed nucleon pole terms. Indeed, assume we have 
the diagrams (12) and (14) of Fig 1 substituting the left A by a nucleon. For small 
momenta of the outgoing pion, on the left of the diagram we would have for the sum 
of the two nucleon propagators 


G [a) + G {b) ~---+---~ 0 . 

M + m n — M E(p') — m, r — E(p) 

This cancellation would not occur if we had a resonance in the intermediate state since 
there is a difference of masses and now we would have 



cy+cf ~ 


M + 771 *- Mr E(p) - 7n n - Er(p) 


m-n 


AM 


-m-n 


AM 


• (32) 


We have also omitted terms involving the exchange of vector mesons in diagrams of 
the type of (l)-(3) in Fig. 1. We estimate their contribution to be small as follows. Take 
for instance the equivalent of diagram (1) substituting the virtual pions by virtual p 
mesons. A general expression for the ppnii Lagrangian is given in [^4|]. For our purpose 
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we can get the low energy estimate by using the analogy of the minimal coupling 
generation of the 77 coupling to 7T7i, in which case we obtain a Lagrangian of the type 


£ 


f- 

J P 


ppTTTT 2 (fV 




(33) 


with f p = 6.1, the pmr coupling. The diagram (1) with two virtual p mesons is now 
readily evaluated and we find at the end that its strength is about the same as the 
corresponding diagram with two virtual pions. Short range correlations also help to 
further reduce the pp diagram versus the nn one. The fact that these meson-meson 
terms are found to be small altogether suggests that this should also be the case with 
the inclusion of the p meson in the analogous diagrams. 

As for the u exchange, one can think of a diagram like (1) in Fig. 1 substituting 
one virtual pion by an u meson. The coupling of the u to three pions can be found in 
23J and is of the type 


oc (34) 

where Pi are the momenta of the pions and the polarisation vector of the aWe can 
expect, in principle, contributions similar to the p meson case. However, at thresh¬ 
old this term will identically vanish since only the zero components in Eq. (34) will 
contribute and they are contracted by an antisymmetric tensor. Hence we should ex¬ 
pect contributions from this term comparatively smaller than the one from p meson 
exchange discussed above and negligible for the total amplitude. 

5 Results and discussion 

In the first place let us look at the reaction pp —> pp7i + 7i~. We show the cross section as 
a function of the energy in Fig. 3. We have separated the contribution of several blocks 
of diagrams in the figure. They are calculated using set 1 of parameters c*, c^, but the 
total contribution is given for both sets 1 and 11 ( solid lines ) in order to give an idea of 
the theoretical uncertainties. Although the sum of the terms is done coherently, there is 
in fact little interference in the total cross section. The short-dashed curve corresponds 
to chiral terms, diagrams (l)-(3) of Fig. 1 . As we can see, this contribution is negligible 
in this channel. The dash-dotted curve corresponds to the diagrams (9)-(15) involving 
only A excitations. We see that this contribution is much larger than the former one. 
At low energies it gives a negligible contribution to the cross section but it rises steeply 
as a function of the energy and becomes dominant at large energies. Among all these 
terms, AA excitation mechanism of diagram (12) is the largest above T p — 1 GeV. 
The long-dashed curve stands for diagram ( 8 ) exciting N* and A consecutively. We 
can see that this term is more relevant than the set of A terms at low energies. Finally 
we show in the long-short dashed line the contribution of the set of diagrams involving 
one N* excitation followed by a two-pion decay in s-wave, diagrams (4)-(7). We can 
see in the figure that this gives by far the largest contribution at low energies. The 
sum of all contributions is given by the solid lines. These lines, corresponding to the 
two acceptable sets of clj parameters, differ by about a factor two at low energies 
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and about 30 % at T p ~ 1 GeV. This sets the level of the theoretical uncertainties 
in this reaction. Another reading of these results is that the pp —> pp / K + n~ reaction 
is a more sensitive tool to the N* —> N couplings than the ttN —> tttiN 

reaction and could be used to put stronger constrains on them. In the next section 
we shall see how these results can be improved. In any case these results have to be 
seen with the perspective that by omitting the N* terms the disagreement at energies 
below T p = 900 MeV is larger than two orders of magnitude. The A A mechanism has 
been used in connection with the ABC effect in the np —> d + X reaction [E5| , but at 


energies corresponding to T p > 1200 MeV in the elementary reaction. The strength 
of this term at these energies is sizable but we get two other sources of contribution 
from N* excitation, as we have shown, which have about the same strength as the AA 
term. We should mention, however, that the consideration of short range correlations 
has decreased the contribution of the AA term with respect to the n exchange alone 
by about a factor three. 

In order to show the relevance of the findings of the isoscalar excitation of the Roper 
we show in Fig. 4 the contribution of the N* terms, diagrams (4) and (5) and show 
there three curves. One with the contribution of the term if we assume a correlated 
7T T p exchange in the T=1 channel (V^, VJ terms), short-dashed curve, another one 
with the contribution assuming only “a” exchange , long-dashed curve, and the third 
one the results with the sum of the two exchanges, solid line. As we can see, the 
contribution with “a” exchange is about one order of magnitude bigger than the one 
obtained with the correlated 7r + p exchange. This shows the importance of the novel 
findings on the isoscalar Roper excitation in order to understand the 2n production in 
the pp reaction. This also gives hopes that the permanent problems in the ABC effect, 
tied to the poor angular dependence provided by the AA model 
solution in the light of this new interpretation of the pp 

m 


261 could find a 


In Ref. 


ppir + n reaction. 

a different approach was followed based on the one pion exchange 


model and two mechanisms: one with two pions produced from the same baryon line 
and another one with a pion produced in each baryon line. The ingredients needed 
there, the tiN —> ttttN and nN —> nN amplitudes, were taken from experimental 
cross sections, making several assumptions on how to extrapolate them off shell and 
summing incoherently the contribution of the two mechanisms. The model was used 
at higher energies that those explored here. Even if phenomenologically one would be 
considering in [^7] the terms discussed here with explicit models, it would only account 
for the 7r exchange in the terms with N* excitation followed by N —> N (Trn^Z^ave or 
N* —> Att, while we have shown that the “cr” exchange is the dominant piece in the 
NN —> NN* transition. Also, we have seen that the indirect effect of the short range 
repulsive NN force weakens the 7r exchange contribution in the 2tc production process. 
Thus, the experience gathered through the years on the NN interaction and the pion 
nucleon and nuclear interaction, together with the recent findings on isoscalar Roper 
excitation, have made it possible the detailed model of the present work, clarifying and 
improving the ideas contained in Ref. [[R| . 

In Fig. 5 we show the results for the pn —> p7nr + ir~ reaction with the same meaning 
as in Fig. 3 and similar results, although with a larger discrepancy than in the previous 
case. 
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In Fig. 6 we show the results off the pp —*■ pmr + n° channel. This reaction is inter¬ 
esting because the N* excitation with N* —> N( y 7nr) r g=^ ljave decay shown in diagrams 
(4) and (5), which were dominant in the pp —> ppn + Ti~ reaction, do not exist now. 
Diagrams (6) and (7) still contribute, but they are very small because they involve 
one N*Ntt p-wave coupling, which vanishes at threshold and also a exchange is not 
allowed. Indeed the 7r + 7r° system can only be in T = 1,2 but not in T = 0. Hence, 
the mechanism that was dominant in the pp —> ppTi + Ti~ reaction at low energies is not 
present here. In spite of that, the agreement with the data is of the same quality as the 
one found before for the pp —> ppTT + 7i^ reaction. Now the dominant terms are those 
exciting A’s. 

In Fig. 7 we show results for the pn —> ppTi~n° reaction. The features are qualita¬ 
tively similar to those in the previous channel but the discrepancies are considerably 
bigger. We note that the strength of diagram (8) is now comparatively bigger with 
respect to A excitation terms than in the previous case. 

In Fig. 8 the results for the pp —> 7imr + 7r + reaction are shown. Here, at high 
energies, A terms are still dominant, but below 1 GeV chiral terms dominate the 
amplitudes. This is the only case where these terms are relatively important. This 
also means that we should accept large uncertainties at low energies since, as we saw, 
there are other terms mediated by p meson exchange which are of the same order of 
magnitude and were neglected here. 

In Fig. 9 we show cross sections for the pp —> nn7i°7T 0 channel. This is again 
a channel where the diagrams (4) and (5) are dominant at low energies, like in the 
pp —> pp7r + 7T^ case. Chiral terms are not drawn since they are below the scale of the 
figure. In this case we overestimate the experimental results by about a factor 2-3 
although the quality of the data is not as good as in former cases. 

With these results we exhaust the experimental data and the isospin independent 
channels. Indeed, using the isospin analysis of Ref. [14|] one can see that there are only 
six independent cross sections. For instance one can deduce an interesting relationship 
with the cross sections which is the following 


2cr(pp7r + 7r ) — a(pnn + n 0 ) — 4a(pp7r°7r°) + 2cr(rm7r + 7r + ) 
+2a(p , mT + n~) — 2cr( y ppTT~7i 0 ) — 4cr(pn7r 0 7r°) = 0. 


(35) 


The NNTnr labels stand for the outgoing particles in the given channel. Incoming 
ones are fixed by charge conservation. We have also calculated the pn —> pmr°n 0 cross 
section with our model, as shown in Fig. 10, for which there are no experimental data 
available. We have checked Eq. (35) independently for the different mechanisms of the 
model as a test of consistency which has been passed successfully. 

The model presented here contains only tree level diagrams. Unitarity is not strictly 
fulfilled. Actually, imposing unitarity with four particles in the final state is less than 
trivial, as evidenced by the enormous difficulties in the case of three body final state 
28| . However, one should bare in mind that as far as we have dominance of a resonant 


term, the important aspects of unitarity are included if the proper resonance width 
is used in the resonance propagator, as we do. Partial unitarization is accomplished 
by the introduction of loops, as done for instance in the ttN —> tittN reaction for the 


12 









chiral terms in Ref. 0. However, we saw that this sector plays a minor role in the 
present reaction in view of the dominant contribution of resonant terms. Unitarization 
schemes, as the one of Olsson p9j, have proved to be successful in the two body 
final state when a resonant term is dominant and there is a small background. One 
multiplies the resonant term by a phase e 1 ^, with 0 small in principle, and demand 
that the resulting amplitude satisfies Watson’s theorem, with the global phase of the 
amplitude in a particular channel equal to the one of the final state. The angles 0 
needed in the problem of jN —> nN are of the order of 10° In order to have 

a feeling for what could be the effects of imposing unitarity in our model, we took 
the pp —> pp7r + 7r _ channel and multiplied the dominant N* term by e*A We see that 
for values of 0 up to 20° the cross section changes at the level of 1%. Rough as the 
procedure appears to be, it gives hints that unitarity is not a thing to worry much 
about in the energy domain studied here. 


6 Final state interaction 


We are going to make a qualitative study of the effect of final state interaction (FSI). 
Since the energy of the incoming particles is large at T p > 800 MeV, we take plane 
waves for the initial state and look at modifications only from the interaction of the 
final particles. Since the low energy region in the pp —» ppir + ir~ reaction is dominated 
by the N*(7nr )sZ^ a „ e contribution, we concentrate on this mechanism alone in order to 
assert the effect of the FSI. For this purpose we substitute 

m = (36) 

by 

r (\ P1-P2+P5+P6 A 

/ d 3 np(r)e\ 2 '/(r) (37) 

where f(r) is the Fourier transform of /(g), 

< 38 > 

The momenta are chosen according to Fig. 11. By taking a monopole form factor 


we find 


F(q) 


A 2 — m 2 
A 2 — g 02 + <f 2 


(39) 



A 2 - ml e _ ir + 


0\2 


~ 9 9 

mt = mt 




0 \ 2 


(40) 


On the other hand p{f) is the pp final wave function, which for low energies can be 
written as 
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<p(r) = e( lkr ) + Jo (k, r) - j 0 (kr) (41) 

where k = ( p 4 — Ps)/ 2 , k = \k\ and Jo (k,r) is the interacting pp relative radial wave 
function with the boundary condition at r —> 00 


jo(k,r) —> e ld °—sin(kr + 5 0 ) 
kr 

which we calculate with the Paris potential |3T]]. Thus we finally substitute 


g02_g*2_ r 


F 2 (q) g 02_g<2_ r 


+4vr / 0 °° r 2 dr j 0 (Qr)\j 0 (k, r) - j 0 (kr))\f(r) 


(42) 


(43) 


and close to threshold 


Q ~ “ = Pi {CM) 
q° ~ E(pi) — M ~ m n 


(44) 


which simplifies much the computations. 

The effect of FS1 is an increase of the cross section at low energies with respect to 
impulse approximation (IA), using plane waves. At high energies, the approximations 
of Eqs. (43), (44) are not good but one should expect that the 1A becomes progressively 
more accurate. However, one must also take into account that some of the effective 
couplings are already chosen in a way that they incorporate FS1 effects. This is certainly 
the case in the correlated tt + p exchange which we have discussed, ft is also the 
case in the NN —> NN* transition dominated by the effective “a” exchange, since 
the empirical coupling was obtained in Ref. JTTj] without explicit inclusion of FS 1 . 


Distortions of proton and pion waves in the (a, a') reaction was considered in an eikonal 
approximation in order to eliminate the a breaking channels. Since the analysis of Ref. 


11| was done at an equivalent T p ~ 1 GeV, we find an approximate way to account for 
FS1 in the N* dominated channels in our present calculation by multiplying the cross 
sections by the factor 


<? FSI {T P ) a IA {T p = 1 GeV) 

a IA {T p ) a FSI {T p = 1 GeV) 1 J 

up to T p = 1 GeV and we do not modify them at energies higher than T p — 1 GeV. 
With these considerations we find the results shown in Fig. 12 for the pp —> ppT\ +, n~ 
reaction , which should be taken as indicative of the role played by FS1. We see that 
the slope of the data is better reproduced when FS1 effects are included. 

7 Conclusions 

We have constructed a model for the NN —> unNN reaction consisting of the terms 
appearing from chiral Lagrangians involving nucleons and pions, plus terms involving 
the excitation of A and N* (1440). In the channels where the two pions can be in a 
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T = 0 state, as 7i + 7i~ and 7r°7r°, we find a dominance of the N* excitation in one nucleon 
decaying into N and nn in T = 0, S-wave. The recent experimental findings about 
isoscalar N* excitation in the (a, a') reaction on proton targets are used here and one 
finds that in the pp —> ppTi + n~, pn —> pnn + n^ reactions the NN —> NN* transition, 
driven by the isoscalar V” exchange, and followed by the N* —> N {nn^Zwave decay 
largely dominates the cross section at low energies. This is an important finding of the 
present work which could not have been asserted prior to the experimental observation 
and analysis of the Roper excitation in the (a, a') reaction. As the energy increases, 
the N* excitation followed by An decay takes also a share of the cross section and so 
does the excitation of a A in each of the nucleons, which becomes dominant at energies 
T p > 1300 MeV. Other terms which are calculated are found to play a minor role. 

A different case is the one of the pp —> prm + 7i 0 channel where the N* excitation 
followed by N* —> N [nn^z^ave forbidden. In this case the successive excitation 
of two A on the same nucleon, the A excitation on each nucleon, the N* excitation 
followed by A7T decay and even the chiral terms (at low energies) share the strength 
of the reaction and one obtains a qualitative agreement with experiment. However, 
these same ingredients used in the pp —> ppn + n~ reaction but omitting the N* —> 
>e decay would give cross sections roughly two orders of magnitude smaller 
than experiment at low energies. This gives us a qualitative idea of the important role 
played by this mechanism in this reaction. 

These new mechanisms for the NN —> NNnn reaction are bound to have repercus¬ 
sions in other reactions. The pp —> pp7r°7r° amplitude with one of the two n° produced 
in one nucleon and absorbed in the other one, gives rise to a box diagram that could 
be relevant for the pp —> ppn 0 reaction close to threshold. Similarly, the isotropic piece 
in pp —> ppir + n- coming from the N* excitation followed by the 2n decay might be 
the clue to a better understanding of the ABC effect, which demands such a highly 
isotropic amplitude in order to interpret the angular dependence p6 
direction should be encouraged. 


Steps in this 
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Appendix 

AMPLITUDES FOR THE pp -> ppn + ir- CHANNEL. 

In this channel, the total amplitude can be expressed as 


Ad (i ' — Ad r 3 r 4 rj r 2 (P3 1 P‘1 1 Pli P 2 ) — Ad(l <-»• 2) + Ad(3 <->• 4,1 2) — Ad (3 <-»• 4) (46) 

where the hrst term in the sum is given below for all mechanisms included in the 
calculation. The subindex stands for the number of the diagrams in Fig. 1. 


■ 1 / firNN \ ( (<? ~ 9l) r3ri (a • (3gi + q 2 )) 

7*4 7*2 

V l 

— ^2 ^ ^ 9 1 ^2 ^ ^ 9 ~ [ 2 ( 91 ^ 2 ) +4(p 5 p 6 ) + 3^t 2 ]| F 7T (q 1 )F n (q 2 ) (47) 

9i — A 4 0.2 ~ A 4 J 


Ada 




4|(gi° + 2p 5 °)| h r3ri { 2 ^ + 4|(q , 2 0 + 2p 5 °)| h r4r2 

(A +P 5) 2 -u 2 


(48) 


Ad 4+5 


Ad 6+7 


Ad 8 




P5°P6° 





f-nNN f 

n n 


((J • Q , 2)r 3 ri( cr ' 92 )r 4 r 2 / T W/ \ 

- ^2 - IHW 2 j 

q-2 


Vt(q 2 )) 


+Vip(q r 2) (^")r 3 ri ' (^r^] T 9crNN9vNN* ^ 2 
X [Dat* (p 5 + p 6 + P 3 ) + £>iV* (P5 + P6 - Pi)] 

0 (not allowed by isospin symmetry ) 


(49) 

(50) 


• 1 f* 9nN*A ) fnNN f 

%—■ 


9 n p, 


FV^(q 2 )(a / 


[l pL 


(°' g) mri (p_ • & 7*4 7*2 
^2 
<12 


KM - Crte)) 


mr\ 


[cr. 


7*47*2 


+ 9aNN9oNN* 


dmri ^r 4 r 2 t-, 2 / 

2 - jVaVh 

92 - mi 


X {2(p 5 P6)£r 3 m [Da(P5 + P 3 ) + 3D A (p 6 + P 3 )] “ *(<? ' [P5 X p^]) r3m 
X [D A (p 5 +P 3 ) - 3D A (p 6 +p 3 )]} Dn*(P5 +Pg +P 3 ) 


(51) 
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i(-) ^ E —Da{P5+P6+Ps) 

V p J p p 


![(<?• q2)r 4 r 2 


9 I & 


[5 i(q 2 ■ [p 5 X p 6 ])S V3V1 - (P5 ■ Pe)(v ■ <f2)r 


+4(pi ■ 92 ) (o' -pe) V3V1 ~ (pi ' <12)(V-P5) V3VI } Kte) - 1 4te)) 

+ [5i5 r3ri ((T • [^5 X Pe])r 4 r 2 (Pb ' P&)^^)r^r\ ' (&)r 4 r2 

+4(a • p 6 ) 

V3ri (o' ' pi) r4r2 - (o' • pi) V3V1 (o' • Pe) 7*4 r 2 ] 14 ( 92 ) } ^a(P5 +Ps) 


g |P5 <-»• P6j ^a(P6 + 


- 4 ^/^ 2 r i { 2 ^ + ^ 2 ° + «}^2 ro^ , -NAT 
*-1-) --fe+p^-p 2 - {2(P5-fe + P 6 ))^n 


3 I P 


-i(a- [p 5 x (gl + Pe>)]) r*3n }^A(P5 +P3) 


{2^-A|(g 2 ° + 2p 5 0 )} ( 5 r . 
(92 + P 5) 2 - P 2 


(2(P6 ' (92 + Pb))Sr 


-i(a- [p 6 x (q 2 + pi)]) 

r3ri } Da (pe + 


. 4 * (r\ 2 \{ 2X i + +2j>« 0) } ^ 


3 V P 


(92 +Pe) 2 - P 2 


(2(P5 ' (92 +Pe))<5 r 


+*(o'- [P5 X (92 +P6)])r 3 n}^A(Pl “Ps) 

1 { 2 7T “ + 2p5°)}5r 4 r 2 ^ 

3 - fe+P 5 ) 2 -^ -^(Pe • + P=»^. 


+i(a • [pi x ( 9 I +pl)]) r 3 ri }Da( P i 
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M.12 — 


?; 1 1 f_ 

3 9 l n 



{2(P6 • (?2 + Pe))<W 2 - i{& ■ [p6 X + Pe)])r 4 r 2 } 


(ft +Pe ) 2 

X {2(p 5 • (g 2 + Pe))^ - *(<? • [P5 X (g£ + ^)]) r3 n} 04(92 + Pe) - 14(92 + Pe)) 

{(P5 ' Pe) ( 25 r4n (5 r3r2 "l - 5<5 r3ri <5 r4r2 ) “t - (<7 • P,5)r 4 r 2 (^ ' P6)r 3 ri 


+2i(<J • [^5 X P6])r 4 r 2 5r 3 n 2 %(<J • [p 5 X Pg] )r 3 r 4 $r 4 r 2 } 14(92 3“ 

X j Da(P4 + Pg)-C>a(P3 + P 5 ) 


(55) 


Alls — 


? 11// 
3 9 l /x 



{2(pe • (?2 + Pfi))^r 4 r 2 + i(<? • [p^ x (q 2 + p^)])r 4 r 2 } 


(?2 +P1) 2 

X {2(p 5 • (cf 2 + P6))<5r 3 n + *(<? • [P5 X + p^)])r 3 n} (14(92 + Pe) - 14(92 + Pe)) 
A {(P5 ‘ Pg)( 2<5 r4n (5 r3r2 -|- 5(5 r3ri (5 r4r2 ) -I - {(J • p,5)r 4 r 2 (^ ' Pe)r 3 ri 


2i(cr • (pg X Pg])r 4 r 2 ^r 3 ri A 2?((7 ' [P 5 X Pg] )r 3 ri4 4 r 2 } 14(92 ~7 

xD A (p 2 - Pq)D a (pi - p 5 ) 


(56) 


A'f 14 — 


.11/7 





{2(Pg • (<?2 + Pfi))^r 4 r 2 + *(<? ■ [p£ X (<£ + P^)])r 4 r 2 } 


9 9 \p) \{qh+Po) 2 
X {2(p 5 • (92 + Pe))<Wi - i(<? ■ [p£ x (q 2 + p 6 )}) r3ri } (14(92 + Pe) - 14(92 + Pe)) 
3” {(P5 ' Pg) (25j- 4 ri 4 3 r 2 “I - 3(5 r3ri S r4T2 ) ( 0 " • P,5)r 4 r 2 (^ ' P6)r 3 ri 

2i(a ' [p ,5 X P6])r 4 r 2 4 3 ri 2%{(J • [p ,5 X Pg]) r 3 ri4 4 r 2 } 14(92 A 


X£>a(P 2 - Pg)D A (P3 + P 5 ) 


(57) 
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4 


X {2(p 5 • (<?2 + Pe))<Wi + ^ • [P5 X («2 + P6)])r 3 n} (V£(g 2 + P 6 ) “ V£(g 2 + Pe)) 


“b {(P5 ' P6) (25 r4ri (5 r3r2 ~t~ 3(5 r3ri S r ^ r2 ) (a • Ps)r4r2(^ ' P6)r 3 ri 


+2i((X ' [ps X P6])r4r-2^r 3 r-i 2z((X ' [ps X pg] )r 3 rj } br((?2 "h Pg) 


xD A (p 4 + p 6 )D A (p 1 -p 5 ) 


(58) 


In these expressions p i, P 2 (p = (p°,p)) are the momenta of the incoming nucleons 
while p 3 , p 4 are the momenta of the outgoing ones, p 5 , p 6 are 7r“ and 7r + momenta 
respectively and gp 2 ) denote p 3 ( 4 ) — Pi( 2 )- Here r 3 , r 4 , ri, r 2 , m = (1,2) are spin indices 
and a sum in m is understood in Eq. (51) V' L and V ^ are the longitudinal and transverse 
parts of the NN interaction as defined in Eq. (22), F n (q) and F„{q) are pion and sigma 
form factors respectively. Finally Djy*(p) and D A * (p) stand for N* and A propagators, 
given by the following expressions 




(59) 


Here E;(p) denotes the width of the resonance ||TT 
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Figure Captions 


Fig. 1. Complete set of Feynman diagrams of our model. 


Fig. 2. Dominant mechanisms in the analysis of (a, a') performed in Ref. [TT|] . (a) A 
excitation in the projectile, (b) Roper excitation in the target. 


Fig. 3. Total cross section for the pp —> ppn + n channel as a function of the incoming 
proton kinetic energy in the laboratory frame. Solid line, total ( line labelled 1 for set I 
and 2 for set II of cj, c * 2 parameters ) ; long-short dashed line, N* —► N {nn^z^avei long- 
dashed line, N* —> An] dash-dotted line, A excitation mechanisms; short-dashed line, 
non-resonant terms from diagrams (l)-(3). The partial contributions are calculated 
with set I. Experimental data are taken from Refs. Jl2|, [13], [14 . 

Fig. 4. N* —► N{nn) r g_*l Jave contribution to pp —> ppn + n total cross section ( set I 
) separated in T = 0 “a” exchange ( long-dashed line ) and T = 1 correlated n + p 
exchange (short-dashed line ) pieces. The solid line gives their sum. 

Fig. 5. Same as Fig. 3 for pn —> pnn + n 

Fig. 6. Same as Fig. 3 for pp —> pnn + n° 

Fig. 7. Same as Fig. 3 for pn —> ppn n° 

Fig. 8. Same as Fig. 3 for pp —► nnn + n + 

Fig. 9. Same as Fig. 3 for pp ppn°n° 


Fig. 10. Same as Fig. 3 for pn —> pnn^n 


O^rO 


Fig. 11. Dominant mechanism in the pp —> ppn + n reaction at low energies with 
final state interaction included. 


Fig. 12. Total cross section for the pp —> ppn + n channel with ( dotted lines ) and 
without ( solid lines ) final state interaction for both sets I and II. 
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